Inflection point as a manifestation of 
tricritical point on the dynamic phase 
boundary in Ising meanfleld dynamics 



Muktish Acharyya 
Department of Physics, Presidency College 
86/1 College Street, Calcutta-700073, India 
E-mail: muktish@vsnl.net 
and 

Ajanta Bhowal Acharyya 
Department of Physics, Hooghly Mohsin College 
PO-Chinsurah, Dist-Hooghly, India 



Abstract: We studied the dynamical phase transition in kinetic Ising fer- 
romagnets driven by oscillating magnetic field in meanfield approximation. 
The meanfield differential equation was solved by sixth order Runge-Kutta- 
Felberg method. We calculated the transition temperature as a function of 
amplitude and frequency of oscillating field. This was plotted against field 
amplitude taking frequency as a parameter. As frequency increases the phase 
boundary is observed to become inflated. The phase boundary shows an in- 
flection point which separates the nature of the transition. On the dynamic 
phase boundary a tricritical point (TCP) was found, which separates the na- 
ture (continuous/discontinuous) of the dynamic transition across the phase 
boundary. The inflection point is identified as the TCP and hence a simpler 
method of determining the position of TCP was found. TCP was observed 
to shift towards high field for higher frequency. As frequency decreases the 
dynamic phase boundary is observe to shrink. In the zero frequency limit 
this boundary shows a tendency to merge to the temperature variation of 
the coercive field. 

Keywords: Ising model, Meanfield theory, Dynamic transition, Tri- 
critical point 
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I. Introduction: 

The ferromagnetic system, in tlie presence of a time varying external mag- 
netic field, remaining far from statistical equilibrium, became an interesting 
object of research over the last two decades [1]. One interesting nonequilib- 
rium response is the dynamic phase transition. This dynamic phase transi- 
tion is widely studied in model ferromagnetic system in the presence of os- 
cillating magnetic field [1]. Tome and Oliveira [2] first observed a prototype 
of nonequlibrium dynamic transition in the numerical solution of meanfield 
equation of motion for the classical Ising fcrromagnet in the presecce of a mag- 
netic field varying sinusoidally in time. The time averaged (over the complete 
cycle of the oscillating magnetic field) magnetisation plays the role of the dy- 
namic order parameter. They [2] found that this dynamic ordering depends 
on the amplitude of the oscillating magnetic field and the temperature of the 
system. Systems get dynamically ordered for small values of the temperature 
and the amplitude of the field. They [2] have drawn a phase boundary (sepa- 
rating the ordered and disordered phase) in the temperature field amplitude 
plane. More interestingly, they have also reported [2] a tricritical point on 
the phase boundary, which separates the nature (continuous/discontinuous) 
of the dynamic transition across the phase boundary. This tricritical point 
was found just by checking the nature of the transition at all points across 
the phase boundary. The point where the nature of transition changes was 
marked as the tricritical point. No other significance of this tricritical point 
was reported. The frequency dependence of this phase boundary was not 
reported earlier for the dynamic transition in Ising meanfield dynamics. 

In this paper, we studied numerically the dynamic transition in Ising 
meanfield dynamics. Here, we confined our attention to study the frequency 
dependence of the dynamic phase boundary. We studied the tricritical be- 
haviour and found a method of finding the position the tricritical point on the 
dynamic phase boundary. The frequency dependence of the position of the 
tricritical point was studied here. We also studied the static (zero frequency) 
limit of dynamic phase boundary. 

The paper is organised as follows: In the next section the model and the 
method of numerical solution is discussed. Section III contains the numerical 
results and the paper end with summary of the work in section IV. 
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II. Model and numerical solution: 

The time (t) variation of average magnetisation m of Ising ferromagnet 
in the presence of a time varying field, in meanfield approximation, is given 
as [2] 

T— = -m + tanh( (1) 

where, h{t) is the externally applied sinusoidally oscillating magnetic field 
{h{t) — hosm{ujt)) and T is the temperature measured in units of the Boltz- 
mann constant {Kb)- This equation describes the nonequilibrium behaviour 
of instantaneous value of magnetisation m{t) of Ising ferromagnet in mean- 
field approximation. 

We have solved this equation by sixth order Runge-Kutta-Felberg (RKF) 
[3] method to get the instantaneous value of magnetisation m{t) at any fi- 
nite temperature T, ho and uj{= 27r/). This method of solving ordinary 
differential equation ^ = F{t,m{t)), is described briefly as: 

m(t + c/t)=m(t) + (^ + ^ + li^-t + t) 
where 

ki = dt- F{t,m{t)) 

k2 = dt- F{t+^,m+^) 

h^dt-F(t+^,m+^ + ^) 

h — r^i -U IM m -U 1932fci _ 7200^2 i 7296A:3 \ 

fv^ — ai- r 13 ' 2197 2197 2197 ' 

k,^dt-F{t + dt,m + ^-8k2 + » - ^ ) 

k,^dt.F{t+f,m-'^ + 2k2-'-^+'-^-'-^) (2) 

The time interval dt was measured in units of r (the time taken to flip a 
single spin). Actually, we have used dt = 0.01 (setting r=1.0). The local 
error involved in the sixth order RKF method is of the order of ((it)^(= 10~^^). 
We started with initial condition m(t = 0) = 1.0. 

III. Results: 

The dynamic order parameter Q{= — §m(t)dt) is time average magneti- 
sation over a full cycle of the oscillating magnetic field. This was calculated 
after discarding the values of Q for few initial (transient [4]) cycles of the 
oscillating field. Finally, the dynamic order parameter Q is calculated as a 
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function of T, ho and /. Now, depending on the values of these parame- 
ters the system gets dynamically ordered {Q ^ 0) or disordered {Q = 0). 
This shows a dynamical phase transition which is a nonequilibrium phase 
transition. We have studied the transition and determined the transition 
temperature Tii{hQ,f ) in a very simple way. For a fixed set of values of Hq 
and / the temperature T is varied (in step AT = 10~^) and Q is measured 
as a function of T. Then we calculated the derivative ^ numerically (using 
three point central difference formula; where the error 0{dT^) [3]). The tem- 
perature, at which the derivative ^ is sharply minimum, is considered here 
as the transition temperature T^. In this way, we obtained the dynamic tran- 
sition temperature Td{ho,f ) for all values of ho. Here, we have changed ho 
(with interval Aho = 0.02) and obtained the dynamic transition temperature 
TdihoJ). 

Now, for a particular frequency /, the plot of T^iho, /) against Hq gives 
the dynamic phase boundary. This dynamic phase boundary separates the 
regions of Q 7^ and (5 = 0. For fixed frequency, it was observed that the 
dynamic transition occurs at higher temperature for lower values of applied 
field amplitude Hq and vice versa. Fig.l shows such a variation. For / = 0.2 
and Hq — 0.5, the temperature variations of Q and ^ are plotted in Fig. 1(A) 
and Fig. 1(B) respectively. From the sharp minimum of ^ (in Fig. 1(B)) the 
transition temperature Td{ho,f ) was found equal to 0.725. The same plots 
are shown in Fig. 1(C) and (D) for ho = 0.3 (keeping frequency / = 0.2 fixed). 
Here, the transition temperature was found to be equal to 0.919. It is clear 
from the figure that the transition occurs at lower temperature for higher 
value of the field amplitude. In this way, the entire phase boundary (i.e., 
Td{ho, f) as a function of ho for fixed / = 0.2) was obtained. 

The dynamic phase boundary was obtained for different frequency /. It 
was observed that for a fixed value of the field amplitude ho that transi- 
tion occurs at higher temperature for higher frequency. This observation 
was shown in Fig.2, for fixed ho = 0.4. The dynamic transition occurs 
at TdihoJ) = 0.845 for / = 0.2 (sec Fig.2(A) and (B)) and it becomes 
Td{ho,f) = 0.908 for / = 1.0 (see Fig.2(C) and (D)). From the figures it is 
clear that the transition occurs at higher temperature for higher frequency. 
We have reported the results of dynamic phase boundary for frequencies / 
= 0.01, 0.02, 0.05, 0.1, 0.2, 0.5 and 1.0. 

In the limit / — > we approach the equilibrium behaviour. In equilibrium. 
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^ = 0. Equation (1) takes the form m — tanh ^™^ ' = 0. This equation was 
solved by Newton-Raphson iterative method of finding the root to get the 
value of equilibrium magnetisation m{h, T). At any fixed temperature T, by 
changing the value of the external magnetic field /i, we calculated the coercive 
field hf. (i.e, the field for which the magnetisation just changes its sign). The 
value of the coercive field was found to depend on the temperature T, i.e., 
the coercive field is a function of the temperature. 

The dynamic transition temperature is also plotted against the ampli- 
tude of the externally applied sinusoidal magnetic field, in the same figure. 
For a fixed value of frequency, the transition temperature Td decreases as the 
value of the field amplitude increases. This gives the dynamic phase bound- 
ary, below which we observed the dynamically ordered {Q ^ 0) phase and 
above which the phase is dynamically disordered {Q — Q). Different dynamic 
phase boundary was obtained for different values of frequency and plotted in 
Fig. 3. It is observed that for same value of the field amplitude, the transition 
temperature Td increases as the frequency increases. So, the phase boundary 
gets inflated as the frequency increases. 

As the frequency decreases (/ — > 0) the phase boundary shrinks and 
ultimately it approches the curve of temperature variation of the coercive 
field (continuous line in Fig. 3). In this context, one may think that the 
temperature variation of coercive field acts as the static limit of dynamic 
phase boundary. We studied the nature of the dynamic phase boundary. 
This dynamic phase boundary changes its curvature from one side to other, 
as one changes the field from lower to higher value. This means, the boundary 
has an infiection point (where the curvature changes from one side of the 
curve to its other side) which was detected by calculating the derivative 
The derivative ^ plotted against ho (for fixed / = 1.0) shows (in 
Fig. 4) a very sharp minimum (at hm) which indicates the inflection point of 
Td — ho curve. This minimum occurs at ho — 0.72. This minimum or the 
inflection point has a great signiflcance. If the value of the fleld amplitude 
is less than h„i, the transition is a continuous one. A typical transition, 
for ho = 0.70(< hm), was shown in the left inset of Fig. 4. On the other 
hand, if the value of the field amplitude exceeds km, the transition becomes 
a discontinuous one. The right inset of Fig.4, shows such a typical transition 
for ho = 0.74(> hm)- The observation shows that the inflection point on 
the dynamic phase boundary acts as tricritical point. To search the location 
of the tricritical point on the dynamic phase boundary, the nature of the 
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dynamic transition has to be studied at several points. The TCP is the 
point, where the nature (continuous/discontinuous) of the transition changes 
from one side to other. But in the method of finding the inflection point on 
the dynamic phase boundary, one can get the exact location of tricritical 
point on the phase boundary, quite easily at least in this case. Tome and 
Oliveira [2] studied this dynamic transition in kinetic Ising model in mcanficld 
approximation and observed the existence of a tricritical point on the phase 
boundary. But they did not report any method to find the TCP directly from 
the phase boundary. Here, we found a method of getting the TCP directly 
from the phase boundary. 

We have also studied the change in position of the TCP on the phase 
boundary as one changes the frequency. Following the same method, by com- 
puting the derivative ^ and plotting it against ho, for different frequency, 
we obtained the position of the minima of the derivatives (for different fre- 
quencies) and hence the position of TCP's. This was shown in Fig. 5 for few 
frequencies. The position of TCP's for different frequencies are shown (by 
big black dot) in Fig.3. It was observed that the position of TCP shifts to- 
wards higher field amplitudes (consequently lower temperatures) for higher 
frequencies. 

IV. Summary: 

In this paper, we have reported our numerical results of the study of the 
dynamic phase transition in kinetic Ising model driven by oscillating mag- 
netic field, in the meanfield approximation. The dynamic phase boundary 
was drawn in the temperature-field amplitude plane for different frequen- 
cies of the applied oscillating magnetic field. The dynamic phase boundary 
was observed to get inflated as the frequency increases. As the frequency 
decreases it shrinks and in zero frequency limit it seems to merge to the 
temperature variation of the coercive field (equilibrium case) . 

The important thing that we observed here is: The tricritical point on 
the dynamic phase boundary is the point of inflection of the phase boundary. 
This observation made the task, of finding the position of TCP on the phase 
boundary, much simpler. We have also observed that the position of TCP 
shifts towards higher field amplitude for higher frequency. 
Acknowledgments: The library facility provided by Saha Institute of Nu- 
clear Physics, Calcutta, is gratefully acknowledged. 
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Fig.l. The temperature variation of Q and for different values of ho but 
for fixed / = 0.2. (A) and (B) for ho = 0.5 and (C) and (D) for ho = 0.3. 
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Fig. 3. The dynamic transition temperature Td{ho,f) is plotted against the 
amplitude of oscillating magnetic field ho taking frequency / as parameter. 
Different symbols represent different frequencies. / = l.O(O), / = 0.5(+), 
/ = 0.2(n), / = O.l(x), / = 0.05(A), / = 0.02(*) and / = O.Ol(o). The 
bullets represent the TCP's on the phase boundary. The continuous line 
represent the variation of coercive field with temperature. 
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Fig. 4. The derivative ^ is plotted against the field amplitude ho for a 
particular frequency / = 1.0. Left inset shows the temperature variation of 
the dynamic order parameter for ho — 0.70 and the right inset shows that 
for ho = 0.74. 
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Fig. 5. The derivative ^ is plotted against the field amplitude Hq for dif- 
ferent frequencies. Different symbols represent different frequencies. / = 
l.O(O), / = 0.5(+), / = 0.2(n), / = O.l(x) and / = 0.05(A). The continu- 
ous line in each case just connects the data points. 
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